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ABSTRACT 

One of the primary means of determining the spin a of an astrophysical black hole is by actually 
measuring the inner radius rin of a surrounding accretion disk and using that to infer a. By comparing 
a number of different estimates of rin from simulations of tilted accretion disks with differing black- hole 
spins, we show that such a procedure can give quite wrong answers. Over the range < a/M < 0.9, 
we find that, for moderately thick disks {H/r ~ 0.2) with modest tilt (15°), is nearly independent 
of spin. This result is likely dependent on tilt, such that for larger tilts, it may even be that rin 
would increase with increasing spin. In the opposite limit, we confirm through numerical simulations 
of untilted disks that, in the limit of zero tilt, r-m recovers approximately the expected dependence on 
a. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — X-rays: stars 
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1. INTRODUCTION 

Classical, astrophysical black holes have only two defin- 
ing properties: mass and angular momentum (or spin). 
As in other astrophysical contexts, the mass of a black 
hole can straightforwardly be determined by the appli- 
cation of Kepler's Third Law, provided the black hole is 
orbited by another observable object and the distance to 
the system is reasonably well known (at least to within 
a factor of sini, where i is the inclination of the binary 
orbit relative to the observer's line-of-sight). 

The spin a = Jc/GM, on the other hand, is a very 
different matter. So far, no direct measure of spin has 
been proposed, leaving only indirect means of inferring 
it. A commonly used method is based on making some 
measure of the "inner radius" of the accretion disk. In 
general relativity, stable circular orbits are not permitted 
all the way down to the "surface" of the black hole, but 
are instead restricted to lie outside the margina l ly sta - 
ble limit rms . Beginning with iNovikov fc Thorn^ (|1973f ) , 
it has commonly been assumed that the accretion disk 
will observe a similar restriction and truncate at • Be- 
cause lias a well-kn own monotonic depen dence on the 
spin of the black hole (|Bardeen et al.|[T972f ). its associa- 
tion with an observed inner radius of the accretion disk is 
then used to infer the spin of the black hole. The observa- 
tion of the inner radius can come either from modeling of 
the c ontinuum emission in the therm ally dominant state 
(e.g. iShafee eTal. 200S : iDavis et all [2006 1 or from mea- 
suring relativistica lly broadened reflection features (e.g. 
IWilms et al.l[200l . 

However, it has been shown that there are many dif- 
ficulties with such procedures. First, one has many 
choices for how to define the effective inner radius r-m 
of the accretion disk, and some of th ese can vary con- 
sider ably (factors of 2-3) from (|Krolik fc HawlevI 
|200^ . Even if one restricts oneself to such obser- 
vationally relevant measures as the "radiation edge" 
(the innermost radius from which significant luminos- 
ity emerges) or "reflection edge" (the smallest radius ca- 
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pable of producing significant X-ray reflection features}, 
signif icant deviations may be observed (jKrolik & Hawle^ 
[200l lBeckwith et al.l[2008l ). For instance lAgol fc Krolik 
(2000l) showed that significant emission can emerge from 
inside rms, thus precluding a direct correlation with rjn. 
Similarly Reynolds & Begclman (1997) showed that the 
reflection edge is not necessarily tied to rms- 

Nevertheless, a strong result for untilted disks (disks 
that lie approximately in the symmetry plane of their 
black hole spacetimes and have their angular momenta 
aligned with the spins of their black holes) is that ri„ 
should gradually decrease with increasing black hole spin, 
i.e. as a/M — > 1. Therefore, it might at least be hoped 
that a study of rin for a large ensemble of accreting black 
holes might yield some information about the range and 
distribution of black hole spins represented in the ensem- 
ble. 

However, we show in this Letter that there is a further 
complicating factor. Based on numerical simulations, we 
show that the effective inner radius of tilted disks (disks 
that do not have their angular momenta aligned with 
the spins of their black holes) yield very different results 
for rin than their untilted counterparts for all dynamical 
measures we have considered. A disk could be tilted for 
many reasons. In stellar mass binaries, the orientation 
of the outer disk is fixed by the binary orbit, whereas 
the orientation of the black hole is determined by how 
it became part of the system. If the black hole formed 
from a member of a preexisting binary through a super- 
nova, then the black hole could be tilted if the explosion 
were asymmetric. If the black hole joined the binary 
through multi-body interactions, such as binary capture 
or replacement, then there would have been no preex- 
isting symmetry, so the resulting system would nearly 
always harbor a tilted black hole. This same argument 
can be extended to AGN in which merger events reori- 
ent the central black hole or its fuel supply and result in 
repeated tilted configurations. Tilted disks are subject 
to addition al torques due to differenti al Lense-Thirring 
precession ()Bardeen fc PettersonI [19751 ). 

Our results indicate that for moderately thick {H/r ^ 
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TABLE 1 
Simulation Parameters 



the simulations, all data are time- averaged over the inter- 



val tr, 



7100M to t„ 



7900M. To ensure all de- 
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0.2), tilted (15°) disks, rjn is independent of spin, or, in 
some measures, actually increases with increasing black 
hole spin. These results could have important conse- 
quences for observational efforts to constrain black hole 
spins. 

2. NUMERICAL SIMULATIONS 

The numerical data in t his Letter were taken 
from simulat i ons pr esented in IFragile et al.1 (|2007l ) and 
IFragile et al.1 ()2009f ). as well as two previously unpub- 
lished simulations. All of the simulations us ed the 
Cosmos-h-h GRMHD code CAnninos et The 
simulations presented in [Fragile ct al. (2007) and the 
two previously unpublished simulations used a spherical- 
polar grid in Kerr-Schild coordinates, with an effective 
grid resolution of 128'^, except near the poles which 
were purposefully u nderresolved. The simulations from 
IFragile et al.1 ()2009l ) used either this same spherical-polar 
grid or a "cubed-sphere" grid also in Kerr-Schild coordi- 
nates, with a resolution of 128 x 64 x 64 on each of 6 logi- 
cally connected blocks that are morphed into segments of 
a sphere. As such, all of these simulations u sed what was 
referred to as "high" spatial r esolution in IFragile et al.l 
(|2007f l and IFragile et al.1 ()2009l ). In all cases, the simula- 
tions were initialized with an analyticall y solvable, time- 
steady, axisymmetric gas torus (jPe Villiers fc HawlevI 
l2003i) . threaded with a weak poloidal magnetic field with 
minimum Pgas/Pmag — 10 initially. The magnetorota- 
tional instalDility (MRI) arose naturally from the ini- 
tial conditions. The simulations were all evolved for 
^ 7900M, or 40 orbits at r = lOM in units with 
G = c = 1. Only data from the final 790M of the 
simulation are used in this Letter, ensuring that the disk 
is fully turbulent and transient effects from the initial 
conditions have died down. The principle variables that 
are studied in these simulations are the spin of the black 
hole, varying over the range < a/M < 0.9, and the ini- 
tial tilt between the disk and the black hole, being cither 
or 15°. The different simulations are summarized in 
Table [11 

3. EFFECTIVE INNER RADIUS 

As pointed out in iKrolik fc Hawleyl (j2002l ) , there are 
many possible definitions of Tin- In this Letter we will 
ignore any that rely on knowledge of the radiation field, 
instead sticking with definitions based on dynamical fiuid 
properties that are easily obtained directly from the sim- 
ulations. We explore four such definitions. 

To best represent the quasi-steady states achieved in 



rived quantities are representative of the disk and not the 
low-density background gas, we apply density-weighted 
shell averaging, i.e. we present quantities of the form 

{p)i 



where 




is the average density of a given radial shell i. 



dt (2) 



3.1. Surface Density 

Perhaps the conceptually simplest definition of r^^ is 
where the surface density of the disk drops significantly 
from some relevant peak value, especially if this occurs 
over some relatively narrow radial range. Indeed, we 
generically find that the surface density drops rapidly 
from a loc al peak at a r ound r ~ lOM [see e.g. Fig- 
ure 13 of IFragile et al.l (j2007l )]. For concreteness we 
take E(rin) = Emax/Se, where S(ri) = {p)i{H)i and 
Emax = max[E(ri)] within the interval rBH l£ fi < 20AI . 
The local scale height of the disk is estimated from 



H 



|Vp| 



(3) 



Note that th i s defi nition differs from what we used in 
IFragile et al.l (|2007f ) and gives slightly larger values for 
H, though the resulting profile is quite similar. 

In Figure [1] we plot the resulting estimates of r-m from 
each simulation. For illustration purposes, we also plot 
a line representing the value of as a function of spin. 
The most obvious result of this figure is that, although 
the untilted simulations (circles) follow the trend indi- 
cated by Tnis (though with systematically higher numeri- 
cal values), the tilted simulations [diamonds) do not. In 
fact, instead of rin decreasing with increasing black-hole 
spin, it increases slightly. We further emphasize this dif- 
ference by fitting simple trend lines to the untilted and 
tilted data separately. Clearly the data are following two 
very different trends. This is crucial since it means that 
any estimates of a/M based on the tilted data could be 
very wrong. For instance, the 915H simulation (with a 
spin of a/AI — 0.9 and a tilt of 15°) produces a value for 
Tin that is almost identical to the OH {a/M = 0) simula- 
tion. 

3.2. Inflow Time 

Another possible definition for Tin is where the ra- 
dial infall time ti^{ri) = ri/{V^)i of material moving 
in through the disk becomes shorter than some fac- 
tor of the local orbital time torb('') = 27r/ri, where 
n = (M/r3)i/2/[i + a(M/r3)i/2] ~ ^Kep is the lo- 
cal orbital angular frequency. Specifically, if we take 
iin(nn) = 3iorb(nn), then Figure [5] shows the resulting es- 
timates of Tin. The factor of 3 in our definit ion of fin(nn) 
is rather arbitrary [Krolik fc: HawlevI (|2002f ) used a value 
of 10], and is simply chosen for convenience. Clearly if a 
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Fig. 1. — Plot of the effective inner radius rin of simulated 
untilted (circles) and tilted (diamonds) accretion disks as a 
function of black-hole spin using the surface density measure 
^(^in) = Sm ax/3e. Open symbols use the spherical-polar grid from 
[Fra gile et al. (2007) and filled symbols use the cubed-sphere grid 
from Fragile et al. (200Sj). The solid line is the marginally stable 
orbit Tnis, and the light, dashed lines are illustrative trend lines of 
the untilted and tilted data. 




Fig. 2. — Plot of the effective inner radius rin of black-hole spin 
using the inflow time measure fin(nn) = StorbCnn)- The symbols 
have the same meanings as in Fig. [T] 

smaller coefficient is used, all of the data points in Figure 
[2] would shift to smaller radii, closer to the event hori- 
zon; whereas if a larger coefficient is used, all of the data 
points would move to larger radii. A similar caveat ap- 
plies to all of the normalizations used in this Letter. The 
key point is that the trend lines would remain largely un- 
changed, and it is the trend lines we are most interested 
in, not specific numerical values of Tin- 

Looking at the trend lines in Figure [21 we again find 
with this measure of rjn that the untilted simulations 
closely follow the prediction of Tms- The tilted simula- 
tions again deviate from this trend; instead they now 
show an almost flat trend line. By this measure, for 
tilted disks appears to be independent of spin. Of course, 
this again means that simulation 915H produces a value 
for Tin that is very close to the OH simulation. 

3.3. Turbulence Edge 

IKrolik fc HawlevI (j200 2^ suggested another possible 
definition of rin as the radius where the local turbulent 
velocity T4urb becomes short compared to the local infall 
velocity (V^). Normally in the bulk of the disk, MRI 
driven turbulence produces velocity fluctuations that are 
substantially larger on average than the mean infall ve- 
locity. However, as disk material approaches the black 



Fig. 3. — Plot of the effective inner radius rin of black-hole 
spin using the turbulence edge measure (Vtuih/ {y^))(i'ln) = 
(Vturb/(V"^))max/2e. The symbols have the same meanings as in 
Fig. |T| Note that the values for simulations 515H-S and 515H-C 
nearly overlap. 

hole, its radial infall velocity increases more rapidly than 
its turbulent velocity. For this measure we define the 
inner radius as (Vturb/(^''))(nn) = (Vturb/(V^'~))max/2e, 
where (Vturb/(^^'^))max is defined over the interval vbh < 
Ti < 15M. For our purposes we take 14urb(''i) 
where 



K-ms(''i)i 

(4) 



We see in Figure [3] that this definition of Tin gives sim- 
ilar results to the previous two. Again, ri„ follows the 
expected trend for the untilted simulations, but appears 
to increase with spin for the tilted ones. 

3.4. Magnetic Field Structure 

As a final possible definition of Hn, we consider the 
radius where the local magnetic field structure transi- 
tions from MRLdriven turbulence toward an evolution 
better described as fiux-freezing. Although a purely tur- 
bulent process would not lead to any correlation, the lin- 
ear properties of the MRI and the consistent orbital shear 
of the disk impose some modest cor relations between 
B'' a nd B'^ , such that a^ag 0.001 (jKrolik fc HawlevI 
l200l . where 



*mag 



\u''uf\B\'^ - B''B^\ 



(5) 



is the dimensionless magnetic stress parameter. As the 
flux-freezing limit is approached, correlations between 
the magnetic field components become stronger, and 
Q^mag ^1- In this region, the transfer of energy to 
smaller scales becomes slower than the infall time. For 
our purposes, we take amag(?'in) = 2e(Q!mag)min, where 
(amag)min is defined over the interval rBH < ri < 20M. 
We see in Figure |4] that this definition of Tin gives sim- 
ilar results to the previous ones. Again, rjn follows the 
expected trend for the untilted simulations, but appears 
to be nearly independent of spin for the tilted ones. 

Note in Figure[4]no datum represents the cubed-sphere 
simulation 515H-C. This is because it is very difficult on 
the cubed-sphere grid to properly recover the B^ and 
B'f' components, especially close to the black hole for 
tilted disks. Therefore, we were not able to recover a 
meaningful value for amag much inside 8M for that 
simulation. 
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Fig. 4. — Plot of the effective inner radius rjn of black-hole 
spin using the magnetic field structure measure «mag('"in) = 
2e(amag)min- The symbols have the same meanings as in Fig. 

m 

4. DISCUSSION 

We have explored four possible measures of the effec- 
tive inner radius of simulated black-hole accretion disks, 
based on surface density, inflow and turbulent velocities, 
and magnetic field structure. With all four measures, we 
found two consistent trends: 1) for untilted simulations 
Tin closely follows the analytic form of r^s': whereas 2) for 
15° tilted simulations rin shows a flat or even increasing 
trend as a/M 1. 

The fact that Figures [l]-|4] all show remarkable similar- 
ities suggests that our results are not due to a poor choice 
of definition for Also, although we have ignored the 
dependenc e of on inclinat ion in making our plots. 
Figure 5 of [Fragile et al.l ()2007[ ) shows that the deviation 
would be quite small for i = 15°, as would be appropriate 
for this work. 

We can get some estimate of the uncertainties associ- 
ated with our results by comparing data obtained from 
the two 515H simulations, one done on the spherical- 
polar grid and one on the cubed-sphere. For the first 
three measures of rin, where we were able to get estimates 
using both simulations, they agree remarkably well. Fur- 
thermore, for most of the measures the trend lines fit the 
data well, suggesting small uncertainties. 

We conclude that the discrepancy of rin between tilted 
and untilted simulations is a robust result, and must, 
therefore, be rooted in some physical process. The 
most likely agent of such a process would seem to be 
the standing shocks aligned along the line of nodes be- 
tween the black hole sy mmetry plane and disk midplane 
(|Fragile fc BlaesI I2008D . These features are not present 
in simulations of untilted disks. Such non-axisymmetric 
shocks can be quite efficient at extracting angular mo- 
mentum from the gas flow, which would reduce the effect 
of the black hole spin on rin . 

In this work we have only considered two values of ini- 
tial tilt: Pi = and 15°. It stands to reason, however, 
that as Pi the inner radius would approach the un- 
tilted values obtained in this work. On the other hand, 
for Pi > 15°, it may well be that the discrepancies in 
Tin would be even larger. In such a case, might truly 
increase with increasing spin. 

The results in this Letter then imply that, at least for 
moderately thick accretion disks {H/r > 0.1), measure- 
ments of Tin are not reliable predictors of a, unless one 



can independently confirm that the disk is not tilted. 
This could pose a problem for recent attempts to de- 
termine spin using so-called "Hard" state observations 
(jMiller et all 120061 : iReis et al.l 120081 ) . where the flow is 
generally taken to be thick. 

For very thin disks {H/r < 0.01), the effect of tilt 
is expected to be much different. Such disks are ex- 
pected to be subject to t he Bardeen-Petterson effect 
(jBardeen fc Pettersonlll975l ). where the midplane of the 
disk at small radii would be aligned with the symme- 
try plane of the black hole through the competing action 
of Lense-Thirring precession and viscous diffusion, while 
possibly leaving an outer disk that is tilted at large radii. 
Since the inner disk would be aligned with the symmetry 
plane of the black hole, the effective inner radius should 
be the same as for an untilted disk. This would apply 
for disks in the "Soft" or "Thermally-Dominant" state, 
with luminosities well below the Eddington limit, which 
are likely to be thin. Attempts to estimate the black 
hole spin by modeling the disk in this state may be unaf- 
fected by the cautions raised in this Letter. However, for 
sources with disk luminosities near or above the Edding- 
ton limit (e.g. Middleton ct al. 2006), the disk should 
again be thick, and the effect discussed here will apply. 
Even if the disk is thin, we reiterate that it is still im- 
portant to independently fit for the inclination of the 
inner disk, since this may not be the same as the binary 
inclination. 

Other methods of determining black-hole spin, such as 
using quasi-periodic oscillation (QPO) frequencies, may 
not be affected by our results. For instance, QPO models 
based on resonant frequ encies attributable to the black 
hole spacetime itself fe.g. lAbramowicz" fc Kluzniakll200ll) 
would likely not be strongly affected. In fact, the weak 
dependence of on a detailed in this Letter actually 
helps one model of low-frequency QPOs based on the 
precession of a radially extended thick disk, by ensuring 
that the maximum frequency would not e xceed ^10 Hz 
regardless of black hole spin, as observed (jlngram et al.l 
2009). However, to date, no consensus model for QPOs 
has emerged, leaving us with few options for unambigu- 
ously determining black-hole spin. 

To end, we emphasize that it would be wrong to con- 
clude from this paper that all tilted black hole accre- 
tion disks should appear to have rjn ~ 6M. Remember, 
the normalizations used in each of our plots were chosen 
simply to give values for the untilted disks that were rea- 
sonably similar to the values for rms- There is no other 
physical motivation for choosing those normalizations, 
and a different choice would have led to different numer- 
ical values for rjn. Further, we can not really say how the 
dynamically determined in this work would compare 
to values of ri„ measured from continuum modeling or 
iron-line fitting without doing more work. Specifically, 
we would need to generate synthetic disk images from 
the simulation data, and measure a radiation or reflec- 
tion edge. Such work is alr eady underway using the rel- 
ativistic ray tracing code of lDexter fc Agoil ([2009'). Even 
then, there is the problem that the original simulations 
did not properly account for all of the important physi- 
cal processes occurring in the disk, notably heating due 
to dissipation of turbulent energy and radiative cooling. 
That will have to await future simulations. 

What can be concluded from this paper is that all tilted 
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disks (with a tilt around 15°) should appear to have the 
same rin, independent of spin. Our results may further 
imply that measurement of a small value for rjn (with- 
out defining small) would require a rapidly rotating black 
hole and an untilted disk, whereas a large value of fin 
would not necessarily mean that the black hole is spin- 
ning slowly; it could simply be tilted. 
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